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We study an optomechanical system in which the mechanical motion of a single trapped ion is 
coupled to a cavity field for the realization of a strongly quantum correlated two-mode system. We 
show that for large pump intensities the steady state photon number exhibits bistable behaviour. 
We further analyze the occurrence of normal mode splitting (NMS) due to mixing of the fluctuations 
of the cavity field and the fluctuations of the ion motion which indicates a coherent energy exchange. 
We also find that in the parameter regime where NMS exists, the steady state of the system shows 
continuous variable entanglement. Such a two-mode optomechanical system can be used for the 
realization of continuous variable quantum information interfaces and networks. 
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I. INTRODUCTION 


Cold trapped ions represent elementary quantum systems that can be isolated from the environment. These trapped 
ions can be brought nearly to rest by laser cooling and both their internal electronic degree of freedom and external 
motion can be coupled and manipulated by external lasers. This unique property makes the cold trapped ions ideally 
suited for quantum optical studies under well controlled conditions. 

Cold trapped ions are one of the most promising physical systems for implementing quantum computation [l| 
and quantum simulations 0, search for new physics bewnd the Standard Model ultra-high precision measure¬ 
ments 0 and extremely accurate frequency standards The long coherence times and individual addressing of 
ions makes them suitable for the experimental implementation of quantum gates and quantum computing systems. 
Recent experimental progress with cold trapped ions include simulating a quantum magnet with two-trapped ions 
Q , demonstration of a phonon laser Q and realization of a quantum phase transition of polaritonic excitations using 
two trapped ions 0. The motion of the center-of-mass (CM) of an ion confined in a Paul trap can be regarded as 
a quantum harmonic oscillator. The interaction between the CM and internal degrees of freedom of the ion can be 
coupled by an external laser light in the Jaynes-Cummin gs ( JC) model Q. Based on this JC model, a fundamental 
logic gate has already been demonstrated experimentally [^, together with nonclassical motion of states of a single 
ion (lll| . 

In the present work, we study an optomechanical system comprising of a single trapped ion coupled to a single 
optical cavity mode. We study the system in a highly detuned limit so that we can eliminate the electronic degree of 
freedom of the ion. We show that coupling of the fluctuations of the mechanical motion of the ion and the fluctuation 
of the cavity field leads to the splitting of the normal mode into two modes (Normal Mode Splitting). The normal 
mode splitting (NMS) demonstrates a coherent energy exchange between the mechanical and the optical mode. It 
is shown that motional and field fluctuations are entangled in the parameter regime where NMS is found. Thus, in 
this hybrid system, quantum information exchange takes place between the ion and the cavity field. Earlier studies 
related to phonon-photon interactions with a trapped ion in a cavity was reported by Massoni et. al. [T^ and Orsag 
[l^ . They showed the squeezing transfer from vibrations to a cavity field. More recently. It was shown theoretically, 
entanglement between motional mode of an ion chain and a cavity field exists where mechanical bistabiltty is found 
M- During the past several years, several theoretical papers on ion-cavity QED context have discussed issues related 
to entanglement generation, preparation of nonclassical states, or realization of quantum gates [Tsl - f^ 

Entanglement is one of the important elements of quantum mechanics as it is responsible for correlations between 
observables . It is now rigorously investigated in performing quantum communication and quantum computation 
with a high degree of efficiency [2^ . A useful and detailed review of the theory of entanglement in systems of continuous 
variables was recently given by Eisert and Plenio in [^ . Some of the interesting experimental demonstrations in this 
direction are given by the entanglement between collective spins of atomic ensembles [30l| . and between Josephson- 
junction qubits [3lj |. 

The simplest scheme capable of gene rating stationary optomechanical entanglement was studied with a single 
movable Fabry-Perot cavity mirror 321, or both movable cavity mirrors [3^ . In fact, entangled optomechanical 
systems have exciting application in realizing quantum communication networks, in which the mechanical modes 
play the role of local nodes where quantum information can be stored and retrieved, and optical modes carry this 
information between the nodes. This type of scheme is proposed based on free-space light modes scattered by a 
single reflecting mirror in Refs, . This enables the implementation of continuous variable (CV) quantum 

teleportation [^, quantum telecloning 351, and entanglement swapping [s^- Further, entanglement in the steady 
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state of a system is significant because it is stationary, i.e., it has a virtually infinite lifetime, and hence could be used 
repeatedly. 


II. THE MODEL 

We consider a single two-level trapped ion placed inside an optical cavity. We assume that the trapped ion satisfies 
the Lamb-Dicke limit. The optical cavity supports a single mode quantized radiation field with frequency Wc. We 
assume for simplicity that the motion of the ion is restricted to the x— direction i.e along the standing wave traveling 
in the x— direction. The trapped ion experiences a dipolar transition, formed by the electronic ground state \g > and 
the excited state |e > at the ionic transition frequency uja, coupling strongly with the single mode cavity frequency 
ujc- The cavity is pumped by an external laser at frequency Up. The dynamics of the system comprising of the single 
cavity mode and the ion is described by the Hamiltonian, 


H — Hcav + Hgi + Hicm + Hint, (1) 

where Hcav = —hAca^a — ih{ri*a — r]a^) is the Hamiltonian for the cavity mode in a reference frame rotating with the 
pump frequency ujp. Here and a are the creation and annihilation operators of a cavity photon, = ojp — ujc is the 
detuning between the laser and cavity-mode frequency, while the frequency 77 denotes the pump strength. The term 
Hei = —hAa<J^cf~ is the internal ionic dynamics, with = ujp — uja the detuning of the pump from the electronic 
transition and a~^, a~ are the usual pseudo-spin operators, (j+ = |e >< e|, cr“ = \g >< e\. The term = hojmb^b 
accounts for the vibrational energy of the ion with Wm as the frequency of the center-of-mass (CM) vibration and b{b^) 
denoting the destruction and creation operators of the CM vibrational motion mode. Finally, the dipolar coupling 
between the cavity and the ion reads in the rotating wave approximation (RWA), Hint = hg{x){a'^a + a~a^), where 
g{x) = Here Q is the coupling parameter directly proportional to the ion-radiation interaction strength 

(i.e the electric dipole moment of the ion). Here x = b"^ + b denotes the dimensionless position operator for the CM 
of the ion. Also fj = ky^h/2mujm is the Lamb-Dicke parameter with k as the wavenumber, m is the mass of the ion. 
<j) accounts for the relative position of the CM of the ion to the standing wave. When the ion is centered at the node 
of the standing wave, (j) = 0. 

We will work in the large ionic detuning limit (i.e large A^) so that we can adiabatically eliminate the ionic degree 
of freedom by putting a~ = This finally leads to the optomechanical Hamiltonian, 


Horn = fiAa^a + fvjJmb^b + hgo{b^ + b) + hgoa'a{b^ + b) — ih{g*a — ga^), ( 2 ) 

where, A = (^ - A^), go = 

III. OPTOMECHANICAL BISTABILITY 

One of the hallmark of cavity optomechanical system is the optical bistability via optical pumping, which has been 
observed in various systems experimentally. In the current system, it is the nonlinearity due to the coupling between 
the cavity mode and the CM motion of the ion that is responsible for bistable behavior. In order to demonstrate the 
bistable behavior, we start with the Heisenberg equation of motion for the cavity mode and the mechanical mode of 
the ion, 


d = —iAa — igo{b^ + b)a — 70 — 77 , (3) 


b = —iujmb — igoo)a — igo — Td, (4) 

where 7 and T is the cavity damping and the mechanical damping respectively. Substituting the steady state value 
of the mechanical mode bg into the steady state expression for the cavity mode o^ leads to the cubic equation for Og, 
which shows typical bistable behavior. 


= \V?, 


“b G2(A -|- G)x‘^ g + ( 7 ^ -f (A -|- G)^)Xm,s 


(5) 
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Figure 1: Schematic diagram of a single trapped ion in an optically pumped cavity. The mechanical mode of the ion is coupled 
the optical mode. 

where, Xm,s = l^sP and G = ^ . Clearly it is seen that the bistability in the intracavity photon number 

disappears when we have (j) = mr and (/> = mrj2. 



A/r 


Figure 2: Plot of steady state intra-cavity photon number |as|^ as a function of A/y for cj) = 7r/4 (solid lines) and (j) = 0(dashed 
line). The parameters taken are rj/'y = 150, ?) = 1, F/y = 0.01, A^/y = —10, ojml'y = 10. 

In Fig.2, we represent |asp as a function of A/y for (j) = 7r/4 (solid lines) and (j) = 7r/2(dashed line).The system 
shows bistability only for higher pump intensity and (j) = tt/A. Bistability in the intracavity photon number disappears 
when we have (j) = 0. For pump rates higher than a certain critical value, we find three steady state solutions for 
the mirror displacement, with two of them being stable.The system prepared below resonance will follow the steady 
state branch until reaching the lower turning point, where a non-steady-state dynamics is excited. This dynamics is 
governed by the time scale of the mechanical motion of the ion because the cavity damping is almost two orders of 
magnitude faster (F « k). 


IV. FLUCTUATIONSrNORMAL MODE SPLITTING 

Here we show that the coupling of the fluctuations of the ion motion and the cavity field fluctuations leads to the 
splitting of the normal mode into two modes (Normal Mode Splitting(NMS)). The optomechanical NMS however 
involves driving two parametrically coupled nondegenerate modes out of equilibrium. The NMS does not appear in 
the steady state spectra but rather manifests itself in the fluctuation spectra of the ion displacement. We now derive 
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the Heisenberg-Langevin equations for the canonical variables and introduce the input noise operators ain{t) and 
bin{t). We shift the canonical variables to their steady state values, a —>■ + a and b ^ hg + h and linearize to obtain 

the following Heisenberg-Langevin equations: 


a = (-iA - 7)0 - ig 2 {b^ +b) + 


( 6 ) 


b = {—iuJm — r)a — ig2{a^ -b a) -I- V^bin, 


(7) 


where A = A -b gi, = go{bs + b*g), 32 = ffolosl- We will assume throughout our work that L < 7 . The input 
noise operators for the optical field satisfies: < ai„(t) >= 0 , < aj„(t')ai„(t) >= na5{t' — t), and < ai„(t')aj„(t) >= 
(ua + — t). The noise operator bin{t) is considered to be classical thermal noise input for the mechanical degree 

of freedom of the ion which satisfies: < >= 0, < >=< >= nbS{t' — t). The quantities 

ria and rih are the equilibrium occupation numbers for the optical and mechanical degrees of freedom respectively. 
The equations (6) and (7) and their Hermitian conjugates constitutes a system of four first order coupled equations. 


The stability of the system is given by the Routh-Hurwitz criterion which is 2g2 < 


(A 2 72)u;^/|A| 


LOm » 7 and |A | « LOm- 
= sj is given 


The steady state displacement spectrum Sb{uj) for the mechanical quadrature x/xo 

by, 


UJm for 
--b + b\ 


Sbiuj) 


27r 


HxHl" 2rnfc- 


(A 2 . 


■ 7 ") 


2Au!n 


Tbico) 


( 8 ) 


X ^(w) = -b 2uJm^b{oj) -oj"^ - za;( 2 r -b r{,(w)), 


(9) 


Vlb{uj) = gl 


w — A 


w -b A 


(w — A )2 -1-72 (w -I- A )2 -I- 72 


( 10 ) 


^ f _^ mT 2ujjri-^ \ 

w v(u; — A )2 -I-72 A) 2 -(-72 / 


( 11 ) 


Here x(<^) is the mechanical susceptibility of the ion motion. 

Figure 3 shows the plot of the displacement spectrum Sbito) of the ion versus normalized frequency and the nor¬ 
malized effective detuning (A/ 7 ). We clearly observe the usual normal mode splitting (NMS). The NMS is associated 
with the mixing between the mechanical mode of the ion and the fluctuation of the cavity field around their steady 
states.The origin of the fluctuations of the cavity field is the beat of the pump photons with the photons scattered from 
the trapped ion. The mechanical mode and the optical mode forms a system of two coupled oscillators exchanging 
energy. An important point to note is that in order to observe the NMS, the energy exchange between the two modes 
should take place on a time scale faster than the decoherence of each mode. On the negative detuning side, we did 
not observe any NMS due to the onset of parametric instability. NMS implies cooling of the mechanical mode. 


V. CONTINUOUS VARIABLE ENTANGLEMENT 

In order to show that the steady state of the system shows continuous variable (CV) entanglement between the 
cavity field and the mechanical mode of the trapped ion in an experimentally accessible parameter regime, we begin 
with the equations of motion for the position and momentum quadrature of the optical and the mechanical mode of 
the system. 


Aa = APa - jXa + 


( 12 ) 
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Figure 4: Plot of logarithmic negativity En of the coupled ion-cavity system versus the normalized detuning (A/ 7 ) for rih = 10. 
We have taken F = O.OI 7 , ujm = IO 7 , go = — 3 .O 7 . 


Pa = -KXa - iPa " + y^P^, 


(13) 


Xb = ujmPb - rXb + V^xir, 


(14) 


Pb = -UJmXb - TPb - 2g2Xa + V^P^, 


(15) 


where, Xa = , 
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The system of above linearized equations of motion can be written in the following compact form, 


R{t) = MR{t) + N{t), (16) 

where, R{t) = {X},, Pf,, Xa, Pa)'^ and the noise vector N{t) = Here M is the 

drift matrix, given as: 


( r u)'YYi 

0 

0 

\ 

cjjTi r 

252 

0 


0 0 

-7 

A 


\-252 0 

-A 

-7 

) 


(17) 


The formal solution of Eqn.(17) is given by R{t) = F{t)R{0) + J* dsF{s)N{t — s) with F{t) = The system 
reaches a steady state only if it is stable, which is possible when all the eigenvalues of the drift matrix M have negative 
real parts so that P(oo) = 0. Consequently, one gets the following equation for the steady state correlation matrix 
(CM): 


MV + VM^ = -D, (18) 

where, V = dsF{s)DF{s)'^ and D is the diffusion matrix given as D — Diag\r{2nb + l),r(2n6 + 1),7,7]. In 
order to study the conditions under which the optical mode and the vibrational mode of the ion are entangled, we 
compute the logarithmic negativity Em which is for the CV case is 


Piv = waa:[0. In 2/r ], (19) 

where, = ^[^~ ~ with A = J^(V) and = det{X) + detiY) — 2det{Z). Here we have 

used the 2x2 block form of the CM as: 


r = (20) 

A Gaussian state is entangled only if <1/2 or Adet{V) < '^{V) — 1/4 and it is equivalent to Simon’s necessary 
and sufficient entanglement non-positive partial transpose criterion of the Gaussian states. 

Fig.4 shows the plot of logarithmic negativity Em of the coupled ion-cavity system versus the normalized detuning 
(A/7). The stationary entanglement between a single driven optical cavity field mode and a mechanical mode of the 
ion via radiation pressure is clearly seen in the same parameter range where the NMS exists. Hence, entanglement and 
mechanical cooling exists together. A similar observation was also noted for a zigzag ion chain .This observation 
indicates a strong correlation between the two modes when a coherent energy exchange is taking place. The stationary 
entanglement between a s ingj le driven optical cavity field mode and a mechanical resonator via radiation pressure has 
already been analyzed in (32l |. 

Now, we discuss the optomechanical quantum nondemolition (QND) phonon and photon detection scheme. For 
the optomechanical QND phonon detection, the cavity mode is pumped with a laser and the transmitted signal is 
measured using a photodetector. The transmitted signal from the cavity mode should be filtered first before measuring 
using a photodetector. Moreover, in order to detect the photon number within its lifetime, it is also required that the 
measurement time should always be less than the inverse of the product of cavity decay rate and the photon number. 
Here, the photon number in the detection mode corresponding to different phonon states is to be studied as described 
in |38l| . This detection mode photon number should follow the time evolution of the mechanical mode. In this way, 
the continuous monitoring of photon counts at the photodetector gives the QND measurement of the phonon number. 


VI. CONCLUSIONS 

In conclusion, we have shown that an optomechanical system formed by a cavity mode and the mechanical mode of an 
ion can coherently exchange energy and produce stationary entanglements in an experimentally accessible parameter 
regime. The steady state photon number in the presence of large laser pump intensity shows bistable behaviour and 
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exhibits a dependence on the relative position of the CM of the ion and the standing wave. Bitability is absent when 
the ion is centered at the anti-node of the standing wave. The coupling of the mechanical oscillator, the cavity field 
fluctuations leads to the splitting of the normal mode into two modes (Normal Mode Splitting) indicating energy 
exchange between the two modes. Such a strongly coupled two-mode system showing steady-state entanglement 
can be exploited for the realization of quantum memories and quantum interfaces within quantum-communication 
networks. 
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